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We derive a boson Hamiltonian from a Nuclear Hamiltonian whose potential is expanded in 
pairing multipoles and determine the fermion-boson mapping of operators. We use a new method of 
bosonization based on the evaluation of the partition function restricted to the bosonic composites 
of interest. By rewriting the partition function so obtained in functional form we get the euclidean 
action of the composite bosons from which we can derive the Hamiltonian. Such a procedure respects 
all the fermion symmetries. 
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(3 ■ I. INTRODUCTION 

o : 

' The IBM of Arima and lachello [J, is most successful in describing the low energy nuclear excitations. The bosons 
O of this model are understood as virtual pairs of nucleons, analogous to the Cooper pairs of superconductivity pj- But 
2^ no general procedure to reformulate the nuclear theory in terms of the effective bosonic degrees of freedom has been 
found. 

^ ■ The first attempt in this direction has been performed, as far as we know, by Beliaev and Zelevinsky 3]. But 
^ [ this work makes use of the Bogoliubov transformation which violates nucleon number conservation. Moreover the 
bosonization is achieved only within a perturbation scheme. 
' The first work which relates the IBM to a nucleon Hamiltonian is due to Otsuka, Arima and lachello 4]. These 
^ . authors got exact results for the pairing interaction in a single j-shell. Their result was somewhat generalized 5|, but 
a full solution of the problem has not yet been achieved. 

There are several recipes for bosonization , mostly based on the idea of mapping a fermion model space into a 
' boson space. This requires a truncation of the nucleon space whose effect is in general not easy to control. 
T-H [ In order to avoid the limitations of previous works we try a different approach where we do not assume any property 
• of the composites, other than their dominance at low energy. In particular their structure will be determined only 
' at the end of the calculation. The problem of truncation of the nucleon space will then be traded by the problem of 
. decoupling some bosons from the others, but in a setting where one can hopefully have a better control. 
' To implement Boson Dominance we perform a functional evaluation of the partition function restricted to boson 
composites. In this way we get the euclidean action of these composites and their coupling to external fields in closed 
form. All the fermion symmetries, in particular fermion number conservation, are respected. The bosonization is 
^ therefore achieved in the path integral formalism, and all physical quantities can be evaluated by standard methods, 
r* 1 The first step, necessary also in the derivation of the Hamiltonian, is to find the minimum of the action at constant 
.5^ ' fields. Depending on the solution, one has spherical or deformed nuclei. In the latter case rotational excitations 
: appear as Goldstone modes associated to the spontaneous breaking of rotational symmetry. The notion of spontaneous 
' symmetry breaking survives in fact with a precise definition also in finite systems "?] . We want to emphasize that the 
, closed form of the action opens the way to numerical simulations of fermionic systems in terms of bosonic variables, 
avoiding the " sign problem" . 

To compare with the IBM we can either write the path integral of the latter, or derive the Hamiltonian corresponding 
to our action. We will make here the second choice. But to derive the Hamiltonian we must perform an expansion in 
the inverse of the shell degeneracy. 

Bosonization appears in several many-fermion systems and relativistic field theories. The effective bosons fall into 
two categories, depending on their fermion number. The Cooper pairs of the BCS model of superconductivity, of the 
IBM of Nuclear Physics, of the Hubbard model of high Tc superconductivity Q and of color superconductivity in 
QCD have fermion number 2. Similar composite bosons with fermion number zero appear as phonons, spin waves 
and chiral mesons in QCD. They can be included in the present formalism by replacing in the composites one fermion 
operator by an antifermion (hole) one. Indeed the approach we are going to present can be applied, as far as we can 
see without any conceptual difficulty, in all the above cases, as it has been argued in a brief report of the method 

A different approach to bosonization which also avoids any mapping is based on the Hubbard-Stratonovich trans- 
formation. The latter renders quadratic the fermionic interaction by introducing bosonic auxiliary fields which in the 
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end become the physical fields. The typical resulting structure is that of chiral theories [lOj. In such an approach an 
energy scale emerges naturally, and only excitations of lower energy can be described by the auxiliary fields • At 
present the relation with the present approach has not been fully clarified. 

The paper is organized in the following way. In Sec. 2 we define coherent states of composites and their properties. 
In Sec. 3 we derive the path integral for the composites and find the effective bosonic action. We restrict ourselves for 
simplicity to a nuclear interaction given as a sum of pairing multipoles, but more general forces can easily be included 
and will be discussed in future works. The effective action we derive is, apart from the above limitation, general. In 
Sec. 4 we restrict ourselves to a single j-shell with pairing multipoles and in Sec. 5 we determine the corresponding 
Hamiltonian. In Sec. 6 we report our conclusions. 



II. COHERENT STATES OF COMPOSITES 



Composites of fermion number 2 are defined in terms of the fermion creation operators 

= ^ct^tet . i ^ cL, ct„,. (1) 

Z Z — \ / mi,m2 

nil ,^2 

In the above equation the m's represent all the fermion intrinsic quantum numbers and position coordinates and J 
the corresponding labels of the composites. Composites of fermion number zero can be obtained by replacing one of 
the fermion operators by an antifermion one. The structure matrices Bj have dimension 20 independent of J. Their 
form is determined by the fermion interaction as explained in the sequel, but we assume that they will satisfy the 
relations 

tT{B\BK) = 2Sj^K. (2) 

We also assume them to be nonsingular. Then their dimension is twice the index of nilpotency of the composites, 

which is the largest integer v such that (bj^ ^ 0. It is obvious that a necessary condition for a composite to resemble 

a boson, is that its index of nilpotency be large. But this condition in general is not sufficient, and we must require 
also that 

det(17B^B)" - 1. (3) 

A convenient way to get the euclidean path integral from the trace of the transfer matrix is to use coherent states ^J. 
If we are interested in states with n = n + i' bosons for an arbitrary reference number n we introduce the operator 

'P^=^^^W^ I db*db{b\b)-'\b){b\ (4) 
constructed in terms of coherent states of composites 



\b) = \e^v(^bjb\J, 



(5) 



We would like it to be the identity in the fermion subspace of the composites. Let us see its action on composite 
operators. Let us first consider the case where there is only one composite with structure function satisfying the 
equation 

5t5.ll. (6) 

In order to evaluate the matrix clement {bt\bt-i) we introduce between the bra and the ket the identity in the fermion 
Fock space 

1= ydc*dc(c|c)-i|exp(-c*c))(exp(-cct)| (7) 
where the c*, c are Grassmann variables. We thus find 

/ 1 \ ^ 

(8) 



{bt\bt-i) = j dc*dcE{c\c,blM-i)=[l + ^b;bt-i 
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where 



E{c*,c,b*,b) = exp{~c*c+l-b* cBc+hc*B^c* ] . (9) 



Therefore the action of Vn on the composites 



i^.(^^n-(i-^) (10) 



n 



shows that it behaves hke the identity in the neighborhood of the reference state up to an error of order — rl), 
namely the measure {b\b)~^ is essentially uniform with respect to any reference state. 

It is worth while noticing that in the limit of infinite Q we recover exactly the expressions valid for elementary 
bosons, in particular 
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{bt\bt-i) =1^1 + ^Kbt-ij ^ eMKbt-i), n~^^. (11) 

In the general case of many composites the above equations become 

{bt\bt-i) = [det(l + /3*A_i)]^ (12) 

where = ^j{b.i)tB,j. Then using the condition 13 we find again that V approximates the identity with an error of 
order l/fl 

p\{bir-bir) = I {(bir-Mr- +oiim)). (la) 

Identifying the operator V with the identity in the subspace of the composites is the only approximation we will make 
in the derivation of the effective boson action. 

III. COMPOSITES PATH INTEGRAL 

Now we are equipped to realize our program. The first step is the evaluation of the partition function Zc restricted 
to fermionic composites. To this end we divide the inverse temperature in No intervals of spacing t 

and write 

= tr(7'exp(-iJr))^V (15) 
We will restrict ourselves to a Hamiltonian with interactions which can be written as a sum of pairing multipoles 

H = ct/^oc-^5^ictF^cticF^c. (16) 

K 

The single particle term includes the single particle energy with matrix e, any single particle interaction with external 
fields described by the matrix Ai and the chemical potential /i 

ho = e + M-ti. (17) 

Therefore we will be able to solve the problem of fermion-boson mapping by determining the interaction of the 
composite bosons with external fields. We assume for the potential form factors the normalization 

triF^FK) ^2n. (18) 
For the following manipulations we need the Hamiltonian in antinormal form 

H = Ho-ch^c^ -Y,9kI^Pkc^c^fIc^ (19) 

K 
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where the upper script T means "transposed" and 

h = ho-^gj^F]^FK, Ho^hr{h + ho). 



(20) 



K 



Now we must evahiate the matrix element {bt\ exp(— T_ff)|&f_i). To this end we expand to first order in r (which does 
not give any error in the final r — s- limit) and insert the operator V between annihilation and creation operators 



{bt\exp{-TH)\bt^i) = exp(-iJoT)(6i|7' -ch^rVc^ 
xY^gkr^cFKcV^c^Fj^c^lb). 

k 

Using the identity in the fermion Fock space we find 

{bt\eM-rH)\bt-iy^ = J dc*dcE{c* ,c,b*,bt^i) 
X exp(— iJoT" — c*hTc 



(21) 



lexp gKT Fk c^c* F^c*^ 



(22) 



where the function E{c* , c, b* , h) is defined in (jSJ. By means of the Hubbard-Stratonovich transformation we can make 
the exponents quadratic in the Grassmann variables and evaluate the Berezin integral 



(6t| exp(— ri/)|6t_i) = det i? exp(— iJoT") / JJdaJ^daK 

•' K 



X exp(— a* • a) exp < 2^''^^ 



\ K2 J 



T\-l 



(23) 



where R=l + hT. Setting Tt = (T\- + A*/3t-i) and performing the integral over the auxiliary fields a^* , we get 
the final expression of the euclidean action 

1 



trh+—ti[Hl + p;f3t)+lnTt] 
It 



^(6*,6) = r^|ffo 
t 

"iE^if [tr(rt/3;F|)tr(rtFK/3t_i) + 2tr {VtF],FK 



-tr[rt/3;F^,rtFK/3t-i]+] +-tr [r*/?,* (/3t_i ft^i)] 

where [..,..]+ is an anticommutator. This action differs from that of elementary bosons because 
i) the time derivative terms are non canonical. Indeed expanding the logarithms we get 

tr [ln(a + am + Inr, ] = i Ul^ tUt - \Wtf^mf\ + 



(24) 



(25) 



where \l t j = \ {ft+i — ft) ■ The first term is the canonical one, while the others contain the derivative of powers 
of the boson variables. The canonical form of the first term is due to the normalization of Eq.||2Jl of the structure 
functions, otherwise j3t and /3t-i would not have the same coefficient. Note the difference of the noncanonical terms 
with respect to the chiral expansions, where there are powers of derivatives, rather than derivatives of powers. 

ii) the coupling of the chemical potential (which appears in h) is also noncanonical. Indeed expanding Tt we get 
/itr [PlPt-i — PtPt Pt-i + ■•■) , and only the first term is canonical 

iii) the function T becomes singular when the number of bosons is of order fi, as it will become clear in the sequel. 
This reflects the Pauli principle. 

We remind the reader that the only approximation done concerns the operator V . Therefore these are to be regarded 
as true features of compositeness. 
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The bosonization of the system we considered has thus been accomphshed. In particular the fermionic interactions 
with external fields can be expressed in terms of the bosonic terms which involve the matrix M (appearing in h). 
The dynamical problem of the interacting (composite) bosons can be solved within the path integral formalism. 
This includes the new interesting possibility of a numerical simulation of the partition function which could now be 
performed with bosonic variables avoiding the sign problem. 

Part of the dynamical problem is the determination of the structure matrices Bj. This can be done by expressing 
the energies in terms of them and applying a variational principle which gives rise to an eigenvalue equation. 



IV. THE ACTION IN A SINGLE j-SHELL 



In this paper we restrict ourselves to a system of nucleons of in a single j-shell. Then we identify the quantum number 
K with the boson angular momentum, K = {Ik , Mk), so that the form factors of the potential are proportional to 
Clebsh-Gordan coefficients 

iFjM)m,,,n2^V2nC'^^[^^^, ^=J + l- (26) 

In such a case the structure matrices are completely determined by the angular momentum of the composites and 
the normalization conditions Q Bj = Fj. The points i) and ii) following Eq. H24|l are the only difficulties in the 
derivation of the Hamiltonian which could be otherwise read from the action. We can overcome them by performing 
an expansion in inverse powers of fl. We will retain only the first order corrections, which are of order fl'^ , with the 
exception of the coupling with external fields where they are of order fl^^. In this approximation the first difficulty is 
overcome because noncanonical time derivatives are of order 1/0 and the second one because the only noncanonical 
coupling of the chemical potential of order comes from the only term of the chemical potential of order 0, which 
can be shown to be /i ~ — ^ go independent of the number of bosons. 
The resulting action is 



S{b*,b) = E I E ^5^1 [(^* - + 

+ E E^A ^12 ) IM 

(bis 

hhhh IM ) 



(27) 



where all the b* 's and all the 6's are at time t,t — 1 respectively, and 

'^KiK2 = 7^^^' ( PKiFj^^e) - gi^nSKiK2 



WUisU = (-2.9o + E-9^. lnti[(2/. + l)]i/20<' j j L 



2 



{bjsbiJiM = E '^k%hJ4,M4bhMsbhMi- (28) 

Ms, Mi 

Notice the factor 2 in front of the chemical potential due to the fact that the composites have fermion number 2. 



V. THE HAMILTONIAN 

The Hamiltonian is obtainedlllj by omitting the time derivative and chemical potential terms, and replacing the 
variables &*,6 by corresponding creation-annihilation operators a\a, satisfying canonical commutation relations 

H = E '^IiMil2M2a\^Mi^l2M2 + E E 
I1M1I2M2 hhhli IIM 

{Ki2ish «L) ,,,, («/3 «/J/m} • (29) 
It is easy to check that, due to the symmetries of the 9j symbols, it is hermitian. 
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From the interaction with external fields we get the fermion-boson mapping of other operators 

C^Mc ^tl'(^^iA^i-^-^/U/2) «liMi«/2M, 

/1M1/2M2 
all/,M ^ ^ 

X fl/l A/i A/2 A/3 a/4 A/4 • (30) 

We remind the reader that the above Hamiltonian has been derived under the condition n < < il in a single subshell. 
Therefore if we further assume emim2 — 6'5mim2j the single boson energy matrix is diagonal 0Ji^i2 = {'2e — §1-^^1)61-^1^. 
But the bosonic interactions couple all the bosons with angular momenta for which the 9j symbols do not vanish, 
even if the corresponding potentials do vanish. 



VI. SUMMARY 



We have developed a general approach to the problem of bosonization where we introduce fermionic composites 
without any preliminary mapping of the fermion model space into a bosonic one. Restricting the trace in the partition 
function of the system to the composites we get the euclidean action of the effective bosons in closed form. The only 
approximation made concerns the identity operator in the space of the composites. 

It is perhaps worth while to spend a few words about the nature of this approximation. Indeed it might appear that 
two are the approximations involved. The first one is the restriction of the partition function to composites. This is 
the fundamental physical assumption of Boson Dominance. Then we replace the identity in the composite subspace 
by the operator V, which seems a further approximation. But V differs appreciably from the identity only for states 
with many bosons, states which cannot resemble elementary bosons because of the Pauli principle. We therefore deem 
that the two approximations are essentially one and the same. 

The nuclear dynamics can be studied by the methods of path integrals, including numerical simulations which now 
are not affected by the sign problem. 

To derive the Hamiltonian of the IBM we must make recourse to an expansion in the inverse of the index of 
nilpotency of the composites. In the present work we restricted ourselves to a nucleon model space of a single j-shell 
and to a number of bosons much smaller than the index of nilpotency. Both limitations can easily be removed. 
Concerning the second one, some care must be exercized to respect particle number conservation, as done for instance 
in 0- The first one requires a parametrization of the structure functions according to 

(Sj,A/)mi,m2 = PJ]l32 ^j'l mi 72^2 " ("^1) 
3132 

Now the energies of the bosons are functions of the parameters p. A variational principle applied to these energies 
generates an equation for these parameters. The solution to this equation can in general be found only numerically, 
but the Hamiltonian and the other operators retain their analytic expressions. 
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